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1. 


INTRODUCTION 


In  engineering  applications  the  domains  of  the  problems  under  consideration 
are  often  unions  and  intersections  of  simple  geometrical  objects  such  as  cylinders, 
balls,  cones,  etc..  The  unions  and  intersections  of  these  simple  objects  yield  edges 
and  vertices.  It  is  well  known  that  the  singularities  of  the  solutions  occur  near  the 
edges  and  vertices.  The  singularities  make  the  computation  for  these  problems  on 
the  domains  with  edges  and  vertices  extremely  inefficient  and  inaccurate.  Hence 
precise  description  of  the  singularity  is  not  only  significant  for  the  regularity  theory 
of  partial  differential  equations  on  nonsmooth  domain,  but  also  extremely  important 
for  the  construction  of  effective  numerical  approximation  methods. 

This  paper  is  the  second  of  a  series  devoted  to  the  analysis  of  regularity  of 
solutions  of  elliptic  equations  on  nonsmooth  domains  in  IR^,  and  it  will  concentrate 
on  the  regularity  in  neighborhoohs  of  edges  of  a  polyhedral  domains.  The  typical 
description  of  the  edge  singularity  is  the  asymptotic  expansion  of  singular  functions 
(see  [7,8,9,11,12,13,17,19,20,22]) 

J  S  T 

■‘W  =  EEE  Cjst{x3)i’jst{0)r°'^'^\lnry  +  uq 

j=l  5=0  <=0 

where  (r,  6,  X3)  are  the  local  cylinderical  coordinates,  Cjst  and  ‘^jst  are  analytic  in 
X3  (except  vertices)  and  in  0,  respectively.  Recently  the  classical  weighted  Sobolev 
spaces  and  with  Kontrat’ev-  and  Maz’ya-type  weights  were  used  to  analysing 
the  regularity  of  high-order  derivatives  of  solutions  (see  [18,21,23]).  As  indicated  in 
previous  paper  [15],  these  approaches  do  not  sufficiently  characterize  the  behaviour 
of  solutions  near  the  edges.  The  solutions  u{x)  in  the  edge-neighborhood  is  analytic 
except  at  the  edge,  and  their  derivative  of  order  k>l  may  grow  rapidly  as  x  tends 
to  the  edge  and  as  k  increases.  The  regularity  results  in  terms  of  the  asymptotic 
expansions  and  the  classical  weighted  Sobolev  spaces  are  unable  to  reflect  these 
natures  of  regularity  in  the  edge-neighborhood.  The  classical  weighted  Sobolev 
spaces  W|  and  with  0  <  <  1  are  suitable  only  for  the  regularity  of  lower-order 

derivatives  of  the  solution,  but  not  for  higher-order  derivatives,  for  instance,  k  >  2 
if  the  elliptic  equation  is  of  the  second  order. 

In  this  paper  we  will  analyze  the  regularity  of  solution  in  the  edge-neighborhoods 
in  the  frame  of  the  weighted  Sobolev  spaces  and  countably  normed  spaces  with 
dynamical  weights.  The  theory  of  these  spaces  on  the  edge-neighborhoods  has  been 
well  established  in  previous  paper  [15].  The  regularity  results  in  terms  of  these  spaces 
will  provide  us  with  the  complete  qualitative  and  quantitative  informations  of  the 
derivatives  of  solution  at  all  orders  and  will  lead  us  to  the  exponential  convergence  of 
the  approximation  by  properly  selected  piecewise  polynomial  spaces  (see  [8,9,14,16]). 

Although  the  regularity  results  for  problems  in  vertex-neighborhoods  of  polygo¬ 
nal  domains  are  similar  to  those  for  problems  in  edge-neighborhoods  of  polyhedral 
domains  (see  [2,3]),  it  is  worth  indicating  that  there  are  differences  on  the  substances 
and  approaches.  We  will  elaborate  the  substantial  differences  in  Section  4. 

The  notations  and  definitions  of  various  spaces  will  be  quoted  in  Section  2  from 
the  previous  paper  [15].  The  Section  3  deals  with  the  existence  and  uniqueness 
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of  weak  solution  of  Poisson  equation  on  polyhedral  domain  with  data  given  in  the 
corresponding  weighted  Sobolev  spaces.  The  main  part  of  the  paper  is  Section  4 
in  which  the  regularity  of  solutions  in  the  edge-neighborhoods  will  be  derived  in 
the  frame  of  the  dynamical  weighted  Sobolev  spaces  and  countably  normed  spaces. 
These  regularity  results  for  Poisson  equation  can  be  generalized  to  linear  elliptic 
equation  and  system  of  equations  without  substantial  difficulties. 

2.  PRELIMINARY 

We  shall  quote  the  notations  and  definition  of  the  spaces  which  were  introduced 
in  Part  I  and  will  be  used  in  this  paper. 

Let  n  be  a  polyhedral  domain  in  IR^  as  shown  in  Fig.  2.1,  and  let  Ti,i  e  J  = 
{1,2,  3,  •••,/}  be  the  faces  (open).  Ay  be  the  edge  which  is  the  intersection  of  P,- 
and  Pj,  and  €  M.  —  {1,2,  ••  •  M}  be  the  vertices  of  Cl.  By  Jm  we  denote  a 

subset  {j  €  I\Ajn  G  Tj.}  of  I  for  m  G  M.  Let  £  =  {ij\i,j  G  /,?,•  (1  fj  =  A,-,},  and  let 
Cm  denote  a  subset  of  £  such  that  Cm  =  {ij  G  C\Am  G  A.j}.  We  denote  by  w.j  the 

interior  angle  between  Pj  and  Tj  for  ij  G  £.  Let  P°  =  (J  P,-  and  P^  =  (J  P,-  where 

ieD 

£>  is  a  subset  of  I  and  A/*  =  I\V.  Further,  let  T>m  ='DT\Im  and  Afm  =  for 

m  G  M. 


Fig.  2.1  Polyhedral  domain  Cl 

For  precise  description  of  the  regularity  of  solutions  of  elliptic  problems  in  poly¬ 
hedral  domains.  We  have  decomposed  in  [15]  the  domain  into  neighborhoods  of 
edges  and  vertices  as  shown  in  Fig.  2.2  and  introduced  the  weighted  Sobolev  spaces 
and  weighted  continous  function  spaces,  and  the  countably  normed  spaces  in  these 
neighborhoods.  The  structures  of  these  space  have  been  fully  studied  in  [15]. 

Assume  that  the  edge  A,j  lies  in  the  0:3— axis  and  Ay  =  {(0, 0,  xa)]  a  <  X3  <  b}. 
Then  a  neighborhood  of  Ay  is  defined  as 

(Ay  )  =  {x  G  fijO  <  r  =  dist{x,  Ay)  <  Ey,  a  -I-  5y  <  X3  <  6  -  ^y} 
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0  <  eij,Sij  <  1  are  such  that  flff  =  <^  for  ^  ^  i,j. 

By  Osmi-^rn)  we  denote  a  neighborhood  of  the  vertex  Am 


Osmi^m)  =  {a:€n  1  0<  p  =  dist{x.  Am)  <  dm} 

Here  Am  is  assumed  the  origin  and  p  =  +  Xj  +  dm  €  (0, 1)  is  selected  such 

that  nf<  =  (f>  for  any  £eCm- 


Fig.  2.2  Neighborhoods  of  edges  and  vertices 
(a)  the  neighborhood  (b)  the  neighborhood 

(c)  the  neighborhood  V5^,<ry Ay);  (d)  the  inner  neighborhood  OsmiAm)- 

Osm{Am)  is  further  decomposed  into  an  inner  neighborhood  of  vertex  and  several 
neighborhoods  of  vertex-edge.  For  ij  €  Cm  we  define  a  neighborhood  of  the  vertex 
Am  and  edge  Ay 


Ay)  =  {x  €  Osm{^m)\0  <  <j><  CTy} 

where  <j>  =  <f>{x)  is  the  angle  between  Ay  and  the  radia  from  Am  to  x.  We  assume 
further  that  Ay  lies  in  the  positive  xa-axis.  Then  sin^  =  yfxl  +  x\l p.  dij  €  (0, 1) 
is  such  that  V5„,5y(Ani,  Ay)f|V5„,<T;t<(ATO,  Afc,f)  =  Am  for  all  ij  €  Cm  and  M  €  Cm, 
ij  ^  M.  The  inner-neighborhood  dsm{Am)  is  defined  as 
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Osmi^rn)  =  Osm{Am)\  U 

ij€:Lm 

For  the  sake  of  simplicity,  we  shall  write  or  ZY(Aij)  Vm,ij  or  V(A„,  Aij),  Om 
OT  d( Am)  instead  of  Ueij,Sij{Aij),Vsm,Sij{Am,  Aij)  and  OsmiAm)- 

By  >  0  integers  we  denote  the  usual  Sobolev  spaces  on  Cl  with  the 

norm 

0<|a|</: 

where  a  =  (01,02,03),  I"!  =  «!  +  02  +  as,  D°u  =  usual  we  write 

H“(a)  =  Hj(n)  =  {u€  H\a)\u  =  0  on  r“},  and 

(semi-norm),  and 

The  weighted  Sobolev  spaces  are  defined  individually  in  the  neighborhood  of 
edges  and  vertices. 

For  X  €  Omi  ^  (0, 1/2)  and  integers  Y  >  0  we  define  the  weight  function. 


ni  (^) 


for  |o|  >  /, 
1  for  |o|  <  i. 


and  weighted  Sobolev  spaces  with  integers  k  > 


{d„)  =  {.  I  iiuiikM  ^  \Ko^4v(6.)  <  -) 

We  next  construct  a  weight  function  in  Vm,ij  =  V  (A„„  A.j),  with  integers  £>0, 
=  i^m,  ^ij),  €  (0, 1/2)  and  Pij  €  (0, 1),  as  follows: 


pf^m+]oi\  ^ 

for  Y  <  oi  -t-  02  <  |o| 

^  ^/5m+|or|— ^ 

for  oi  -H  02  <  Y  <  |o| 

.  1 

for  |o|  <  £. 

Then  we  introduce  the  weighted  Sobolev  spaces  over  Vm,ij  with  integer  k  >  £ 


A  weight  function  in  the  edge-neighborhood  Kij  =U  (A,j)  is  defined  as 

asaY/-  A  _  /  for  Oi  -I-  02  >  £-, 

^'3 ~  I  1  for  oi  -1-  02  <  Y, 

with  an  integer  Y  >  0  and  ^ij  €  (0, 1).  Then  the  weighted  Sobolev  space  i^ij) 
is  given  by 
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[  |al<fc 


It  t2  <  OO  >  . 

L  (Wo  )  I 


Here  /3ij  coincides  with  that  for  the  space  ).  As  usual  \D^u\  =  J2  |-D"u| 

’  |a|=^ 

and 


Let  =  {l3m,Tn  €  M.jiij.ij  €  £)  with  €  (0,1/2),  €  (0,1)  be  a  multi¬ 

index.  Then  the  space  (fi)  denotes  the  set  of  functions  such  that  their  restric¬ 
tions  on  Uij,dm,Vm,ij  and  %  belong  to  (Kij),  (Om)  (Vm,,j)  and 

H*  (fio)?  respectively,  for  all  ij  €  C  and  m  €  M,  and 

(^)  ij^C  meM  ^0m[Om} 

meM  ijeCm  “/3m, 

The  regularity  of  solutions  in  the  edge-neighborhood  Uij  will  be  given  in  terms 
of  countably  normed  spaces  with  weighted  Sobolev  norm 

K,  (««)  =  {«!«€  Hj;'  («„)  for  all  k  >  I,  <  era'] 

and  countably  normed  space  with  the  weighted  C*'— norm 

C^.^lUij)  =  [ueC'^iUij)  I  e  C°(Wij)  for  a  with  |q!|  =  |fc|  >  2, 

||r^'-’‘''l“l“^Z)"(u(x)  —  ^(OjO, a:3))||0o^^  ^  <  Cd°'o:! 


||^«(0,0,X3)||cO(;^  )  <  C4kl  for  k  >  o| 

3  ) 


Hereafter  =  (a  -f  ^ij,6  —  Sij),  d°‘  =  a!  =  ai!Q!2fo3!,  C  >  I,  and  d,-  >  1 

are  independent  of  a. 

The  relations  between  the  spaces  (^jj)  and  C'^..  (Idij)  has  been  discussed  in 
[15]  from  which  we  quote  a  theorem. 

Theorem  2.1  B^..  (Uij)  C  C^..^  (%)  C  {l^ij),s  >  0  arbitrary.  □ 

It  is  convenient  to  use  cylindrical  coordinates  (r,  6,  X3)  with  respect  to  the  edge 
Aij  when  we  analyze  the  regularities  of  solutions  in  the  edge-neighborhood  Uij. 
Hence  we  introduced  weighed  Sobolev  spaces  in  the  cylindrical  coordinates 


=  ii« 


^  =  E  ,  <  OO 


and  the  countably  weighted  Sobolev  spaces 

6 


= {u  €  for  all  t  >  £,  -  '^‘^“’1 

where  =  u^oj  002  j;“3. 

The  following  theorem,  which  gives  us  the  relations  between  spaces  in  Cartesian 
coordinates  and  those  in  the  cylindrical  coordinates,  has  been  proved  in  [15]  . 

Theorem  2.2  For  ^  <  2  the  spaces  {Uij)  and  (Uij)  are  equivalent,  and 
the  space  (Uij)  is  equivalent  to  the  space  (Uij).  ° 

3.  WEAK  SOLUTION  OF  POISSON  EQUATION  IN  POLYHEDRAL 
DOMAIN 

Consider  the  Poisson  equation  in  polyhedral  domain  fl 


(3.1) 


-Au  =  /, 


with  f  £  Lp  (fl)  lr<  and  ^  (fi) ,  ^  —  0)  1 

Lemma  3.1  If  /  €  L/?^  (Om)  ,0  <  <  1/2  and  u  €  (^m)? 


Proof.  By  Schwartz’s  inequality 

|/a„,  /vdxl  <  C  ll/llL,„(a.„)  \\p'^”^4v(6m) 

^C'll/llL^(6„)ll^llL^((Pr„) 

Here  we  used  the  fact  that  j3m  G  (0, 1/2)  and  the  imbedding  result  of  Sobolev  spaces 
(see  [1]).  ° 

Lemma  3.2  If  /  €  L^,^  (Uij) ,  0  <  <  1  and  v  6  (ffij).  then 


Proof.  By  Schwartz’s  inequality 
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By  Lemma  5.1  of  [Part  1] 


r  \vf  dx 


,2(l-A,)-2 


<C  J  +  [vf'^  dx 


which  together  with  (3.4)  yields  (3.3). 

Lemma  3.3  If  /  e  ,  ^m,ii  =  {^m,^ij)  with  G  (0,1/2),  and 

Ai  €  (0, 1),  and  u  €  (Vm.ij),  then 


Proof.  By  Schwartz’s  inequality  we  have 


/  fvdx  SII/llL,,„(v.,„)|k''*“  "I 


L"(v„,o)  ‘ 


Let  ^^  =  II2-  ^rn  and  =  1  -  and  Then  by  Lemma  4.2  of 

[15] 


p  (sin  <f>)  V 


T^v  ^  =  /  ^{p^^hl'^dx 


Note  that 


Then  (3.5)  follows  from  (3.6)-(3.8). 

Combining  Lemma  3. 1-3.3  we  have 
Theorem  3.1  If  /  G  L/j  (fl) ,  then  /  G  (fl)^  ,  and 


(H’(q))' 


^  <^II/IIl«(o)- 


Lemma  3.4  Let  G  G  (Idij).  Then  for  v  G  (Uij) 


Fif] 


G  vdsl  <  C  G  TTi.i ,,,  ,  V  Tji 


m^(wo)  r  iiH’(Wo 


8 


Proof.  Let  lAij  —  Qe  ^  Isij^  —  |(^i5^2)  0  < 

I^..  =  (a  +  Sij,b-  6ij),  and  let  7  =  Q^nri.  By  H^;^.((5£) 
Sobolev  space  over  Qs 


■^x\  +  x\  =  r  <  e  I  and 
we  denote  the  weighted 


Then  for  almost  every  X3  €  G  €  HJ;;  (Q.)  and  n  G  (Q.).  By  Lemma  2.11 

of  [2]  we  have 


(3.10)  \G\  |u|  ds<C  ll^llH'(Qe) 

where  (7  is  a  constant  independent  of  0:3,  integrating  (3.10)  in  xz  over  hij  we  obtain 
(3.9).  ° 

Lemma  3.5  If  G  €  (V„,,i,)  with  /?„,  €  (0,1/2)  and  Ai  €  (0,1),  then  for 

v£H^ 


(3.11) 


Gvds 


<C'1|G1| 


(Vm.tj) 


Proof.  Let  Sfj  =  {{(l),0)\0  <  (t>  <  cr,0  <  0  <  uJij}  and  Is  =  (0, 8m).  Then  Vm.ij  - 
5?-  X  h  .  We  may  assume  that  T,  is  in  the  xi  -  X3  plane.  Let  7  =  dSfj  fl  L,-  = 
{ (^,  0)1 0  <  </>  <  <T,  0  =  0}.  Then  by  Lemma  2.11  of  [2] 


(3.12)  f  |G1  |ulds  <  G  llGlljji.i  ||u|ljji/g^1 

where  [Sfj)  and  (5?)  are  the  Sobolev  and  weighted  Sobolev  spaces  over 
namely 


ul 


and 


Ge 


(l>d(f>d0. 


Multiplying  (3.12)  with  p  and  integrating  it  in  p  over  Is  we  get 


(3.13)  jf  |G1  lu|  pdpd(j)  <G  P  II^IIh;;1,(s?  )  II^IIh'(55) 
by  Schwartz’s  inequality 
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By  Lemma  4.1  of  [15]  we  have 

II«IIh“  (55)'*'’ 

(3.14)  <C  I  \p'^^’^-^\Gf  +  p^‘^”'{sm<l>f^'^i\-G4^  +  \-^Ge\^)\dx 

^  ^  Jvm.ij  [  p  psin(l>  } 


^C-IIGII 


Analogously,  using  the  fact  that  ^  (0,1/2)  and  the  imbedding  theorem  of 
Sobolev  space 


(3.15) 


iI'’IIh'(s5)‘*'' 

<C  [ 


2/5m|vp  I  2(1— prn)  (  _ 


1  1 

—Vd,  +  - ; — rVo 

p  p  Sin 


The  combination  of  (3.13)-(3.15)  leads  to  (3.11). 

Lemma  3.6  Let  G  €  (Om^.  Then  for  u  €  (^»")  ^ 


(3.16)  l/  .  GvdS\<C\\G\\jjr..^ 

JrinOm  “/S 


)  Ii^IIh'(o.„)  • 


Proof.  Let  S  =  S\  ^ijeCm  ^ij  where  S  is  the  intersection  of  fl  and  the  infinite 
polyhedral  which  coincides  with  fl  in  the  vertex  neighborhood  O  {Am)  and  S^-j  were 
defined  in  the  proof  for  the  previous  lemma.  Then  Om  =  S  x  !$„■  The  proof  of 
(3.16)  is  similar  to  that  of  (3.11)  except  that 


instead  of  (3.12),  and 


instead  of  (3.14),  and 


|Gl|u|ds  <  C  II^IIh^(s) 

P^^"*  ll<^liH‘(5)  II^IIh^'^(<9„) 


instead  of  (3.15).  □ 

Lemma  3. 4-3. 6  lead  us  to 

Thoerem  3.2  If  G  G  (fl),  then  for  u  G  (fl) 


<  C*  lie'll ll^llH'(n)- 

*■ 

□ 

We  are  now  ready  to  prove  the  theorem  of  the  existence  and  uniqueness  of 
the  weak  solution  for  the  problem  (3.1)  with  /  and  given  in  the  corresponding 
weighted  Sobolev  spaces. 

Theorem  3.3  Let  be  a  polyhedra  in  IR^,  /  G  L/3(fl),  =  G^\rt  and  G^  G 

=  0, 1  with  G  (0, 1/2)  and  G  (0, 1)  for  all  m  G  A4  and  ij  G  C. 
Then  the  problem  (3.1)  has  a  unique  solution  u  G  (0)  (weak  sense)  such  that 
u  -  G  (fl),  and 


(3.17) 


/ 

JdQ 


Gvds 


(3.18)  ll^'llH^n)  -  ^  ll^'llHr’'"'(«) )  ' 

Proof.  We  may  assume  that  —  0.  The  bilinear  form  on  (O)  x  Hq  (fi)  is 
defined  as 

B{u,v)=  /  Vu  •  Vuda;. 

•/ 

Due  to  Theorem  3.1  and  3.2 

FH  =  L  fvdx  -f  J  g^vds 

defines  a  linear  functional  on  (Q),  and 

By  Lax-Milgram  theorem  there  exists  a  unique  solution  u  G  Hq  (D)  for  the  varia¬ 
tional  equation  of  the  problem  (3.1) 


(3.19) 

and 


B{u,v)  =  F{v),  VuGHJ(D), 


which  is  (3.15)  with  g°  =  0.  For  general  case  that  g^  ^  0  (3.16)  can  be  proven  easily. 

□ 
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Remark  3.1  If  |r°|  =  0  and 

(3.20)  fdx  +  g^ds  =  0, 

then  Theorem  3.3  holds  in  the  quotient  space  module  a  constant.  □ 


4,  REGULARITY  IN  NEIGHBORHOODS  OF  EDGES 


We  shall  make  further  investigation  on  the  regularities  of  the  solution  of  (3.1)  in 
neighborhoods  of  the  edges  in  frame  of  the  weighted  Sobolev  spaces  and  countably 
normed  spaces.  We  concentrate  ourself  on  a  neighborhood  ZY12  of  the  edge  A12.  As 
assumed  in  previous  sections,  A12  lies  in  the  xa-axis  and  Uu  =  {(r,  0,  X3)  |  (r,  9)  G  Qe, 
xz  €  Is}  with  Qs  =  {(j',^)  l0<r<£,  0<^<C(;}  and  Is  =  (—1  +  6, 1  —  <5) .  (r,6,xz) 

are  the  cylindrical  coordinates  with  respect  to  A12.  We  further  assume  that  Fi  C 
r°,  r2  C  r^.  For  sake  of  simplicity  we  shall  write  U  =  U^s  =  W12  =  I^e,s  (A12) ,  Q  = 
Qe  etc..  As  in  Section  2  we  denote  D^u  =  D°‘'u  °3  =  u  aj  *  «3  and  'D°‘u  =  P"  u  “3 

=  u,.aigc.2j.“3  with  a  =  (o', 0:3)  =  (q^i?«2,Q;3)  and  ja]  =  [of'l  +  0:3  =  ai  +  a2  +  <^3-  We 
shall  write  H^j^  {U)  =  H^’^  (W12)  and  {U)  =  B^^^  (^12)5  etc. 

4.1  Regularity  of  high-order  derivatives  with  respect  to  the  direction 
along  the  edges 

Lemma  4.1  Let  T  =  {(xi,X3)  |  xi  G  (0,£),X3  G  Is}  and  G  G  H^^^(S').  Suppose 
that  V  G  H^(r)  and  u  =  0  for  xi  =  £  or  X3  =  ±  (1  —  6).  Then 


(4.1) 


j^GAhvds\^<C\\G\\ 


where  AhV  =  |(n  (xi,  xz  +  h)  -v  (xj,  X3)),  (7  is  a  constant  independent  of  G  and  v. 
Proof.  First  we  extend  G  and  v  into  T  =  (-£,£)  x  h  by  symmetric  manner 
with  respect  to  xi-axis.  The  extended  functions  are  denoted  by  G  and  v.  Then 
G  G  H*^^  (T)  and  u  G  Hj  (T).  Further  we  extend  v  in  whole  plane  by  zero  extension 
outside  T,  and  extend  G  in  the  plane  as  well.  Then  the  H^^^-norm  of  G  and 
norm  of  v  are  preserved,  and  AhV  is  well  defined.  Let  G  and  v  denote  the  Fourier 
transformation  of  G  (G  9)  and  v{G  v)-  The  equivalent  norms  of  G  and  v  in  H^^^(IR^) 
are  defined  as 

II^IIh‘/^ir^)  =  1^1 

and 
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Then 


GAhvdS  =  I J,  GAhvdS  =\j^.  GAhvdS 
r  ~  e~'^  —  1 


(4.2) 


For  \H\  <  1 


e-i^h  _  I 

h 


_ : _ < _ _ <  1 

(1  +  £2  +  ,;2)l/2  (1+^2  ^^2)1/2 


and  for  \h(\  >  1 


_  I 


< 


{1+e+v^f^  ~ 


<2. 


Hence  we  have  by  Schwartz’s  inequality 
(4.3) 

I  f  A  “•  1 

<  2 (4  |6f  (1  (4  (i 


Then  (4.1)  follows  from  (4.2)  and  (4.3)  at  once. 

Select  e'  €  (e,  1)  and  S'  €  (0,  S)  and  let  W  =  Us>^s'  —  ^e',5'  (^12)  C  0.  Then 
U'  DU  =  Uegt  and  we  have  the  following  theorems. 

Theorem  4.1  Let  u  G  (fi)  be  the  weak  solution  of  the  problem  (3.1)  with 
f  eJjp  (O)  and  G^  G  (fi) ,  £  =  0, 1. 

(i)  lifeV  {W)  and  G^  G  {U’) ,  ^  =  0,1,  then  G  (U)  ,  and 


(4.4)  - 


(«)/ 


(ii)  If  43  G  Lp,,  {W)  and  Gi^  G  {U'\  then  u,,  G  (W),  and 
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(4.5) 


s  Co  I  ii/iil,,(j)  +  ii/«,iIlj„(«') + ll®1lH;““"'(ti)  II*^»IIhj7''’“'(m’)}  ■ 


(iii)  If  /  =  fi+  f2  with  /i  €  {IC),  and  f2,x3  €  -  G{  +  Gi  with 

Gi  €  H'-'  (W)  and  G^^^  €  then  ^,3  €  H'  (U)  and 


(4.6)  +,5.i 

<  Co  ( ll/llL,,a,  +  ,5,  l|GllHr-'w  +  "l'L("') 


+,5.,  UI®'1h-'(«-)  +  ]  ■ 


Here  Co  is  a  constant  independent  of  0:3,  and 

(4.7)  Mo  =  max|^,  ,  AS  =  6  -  8' ,  Ae  =  e' -  e. 

Proof.  First  we  assume  that  =  0.  Let  AhU  =  ^  (u  (a;  +  hes)  —  u  (x))  with 

h  e  (0,8)  and  63  =  (0,0, 1) ,  and  let  H^(n)  =  |«  €  H^(0)|«  =  0  for  x  € 

By  the  standard  argument  of  difference  quotient  (see,  e.g.  [10]),  we  have  for  any 

«;  €  (H) 

(4.8)  L  V  (A/ju)  Vwdx  =  (A/iVu)  •  Vwdx 

=  —  f  Vu  •  V  (A-hw)  dx 
J  f2 

by  (3.19) 

=  —  /  /(A_/ito)dx—  f  G^{A-hw)dS 
Jq  ari 

=  -  /  f{A-hw)dx-  f  G^  {A-hw)dS. 

Ju>  Jr2C\dU' 

In  the  case  (i),  /  €  {W)  and  G^  €  [W).  We  have  by  Schwartz’s  inequality 
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(4.9)  \l  f  (A-*u,)  dv\  <  C 

by  Lemma  7.2.3  of  [10] 

By  Lemma  4.1  we  have 

iLnsu'^'  s 

by  the  imbedding  inequality  of  Sobolev  space  (see  [1]) 

£  CllG'llijiiy.j  IIHIh‘(«.) 

<c||g-|Ih',„,,IIv»IIi,>,„v 

Here  we  used  the  inequality 

(4-11)  ^  • 

Combining  (4.8)-(4.10)  we  have 

(4.12)  \l  V{A,u)  ■  V»rfx  |<  C  • 

Let  <^i  (xa)  and  <^2  (»')  be  C°°  cut-off  functions  such  that  0  <  (xa) ,  (/?2  (»')  ^  1, 
and 


/  ll 

for 

k3| 

10> 

for 

(4.13)  <^i(a:a) 


Set  Tj  =  ifi  (xa)  (^2  (j™)  and  w  =  rj^AhU.  Then 


,  h  (r) 


for  r  <  e 
for  r  >  s'. 


V  [Ahu)  •  Vw  =  V  {Ahu)  V  {ri^AhU^  =  |»?V  (A/,t«)|^  -f-  2r)  (Ahu)  Vt/  •  V  (A^u) 
which  together  with  (4.12)  yields 


\rjV  {Ahu)\'^  dx 

(4.14)  +2  II^V  (A,u)||l.^^,j  .  ||(A,u) 

<  c  (||/IIl^^,o  +  IIG'^IIhm^o)  (|I»?v(a,u)||l2(^,j  +  ||(a,u) 
+2  ||(Afc«)  V77||l2(^,j  •  ||»?v  (Afeu)||L2(^,j 
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Let  A{f,G')  =  C  +  IIGMIh'im.))'  Then  we  have 


ll(Vau)  II^V  <  i  ||>/V  +  2  ||(Aan)  . 


A  (/,  G')  IKV  (Ait<)|l^a,„,,  <  2  |A  (/,  G')f  +  i  ||.)V  (Aju)|||^.,„, 


A(/,G')ll(Ak<i)Vi;||j_j|j^,j  <  j|a(/,G*)|  +2ll(A/.«)Vi)||j^2j^ 


Substitution  of  these  inequalities  into  (4.14)  gives 


Note  that  rj  =  1  in  U,  |  V?/|  <  CMq  with  Mq  =  max  (^,  .  Then  by  Lemma 

7.23  of  [10]  we  obtain  (4.4)  for  G°  =  0. 

In  the  case  (ii),  U  E  (7/') , ^  =  0,1-  Then  for  any 


(4.15) 

by  Schwartz’s  inequality 


by  Lemma  5.1  of  [15] 


I  fA-hwdx  =  I  {Ahf) 

JQ  JW 


by  Lemma  7.23  of  [10]  and  (4.11) 


^  G  UxAIl^jw) 


Due  to  Lemma  3.2 


(4.16) 


\ir2C\au'G'  (A_»ti))i<S|  <  l/r^paw'  (AjG')u)<iS 


by  Lemma  7.23  of  [10]  and  (4.11) 


<  C  \\AhG  lljjM 
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Combining  this  with  (4.8),  (4.15)  and  (4.16)  we  obtain 


(4.17)  jf  ^  V  (Ahu)  ■  Vwdx  <  C  (i1/x3|IL;3i2(«')  +  II^^3  1|h^;‘^{W'))  ’ 

Then  the  proof  of  (4.5)  is  the  same  as  that  for  (4.4)  except  that  the  inequality  (4.17) 
is  used  instead  of  (4.12). 

(4.6)  for  =  0  in  the  case  (iii)  is  obtained  by  combining  (4.4)  and  (4.5). 

We  now  prove  (4.4)  and  (4.6)  is  general,  i.e.  G°lri  O-  Let  v  =  u  —  G°.  Then  v 
satisfies 

—An  =  /  +  Ai2G°  +  G°2  =  f  1 
^  ^Iri  =  0, 

where  =  ^  +  |^.  If  /  €  L"  (W) .  G'  e  (»'),  then  /  G  V  (W)  and 

G^  €  {W).  Applying  (4.4)  with  =  0  we  obtain  (4.4)  in  general.  If  G 

hp,,  {W) ,  (U') ,  ^  =  0,1  then  f  =  fi  +  ^  with  /i  =  G°2  € 

{td') ,  /2  =  (/  T  Ai2f7°)  €  L/3j2  {^')  f2,x3  €  L;3j2  {t( ),  G  ,Gj.^  G  )• 

Applying  (4.6)  we  have  (4.5)  in  general.  ^ 

The  regularity  of  higher  derivates  in  X3  is  given  in  the  next  theorem. 

Theorem  4.2  Suppose  that  u  G  (fi)  be  the  weak  solution  of  the  problem  (3.1) 
with  G^  G  (fi)  ,^  =  0, 1  and  /  G  L;3  (fi). 

(Bl)  If  G  V{l{'),Gi^  G  =  0,1  ,  0  <  m  <  fc,  then  w^j+i  G 

(U),  and 

(4-18)  ■  . 


(B2)  If  /.»  e  Lft.  (W) ,  G'„  €  («') ,  <  =  0, 1  ,  0  <  m  <  +  1,  then 

u^k+1  G  (W),  and 

(4.19)  .  .  ,  , 


““^3+'  Hb 


Furthermore,  if  /  G  {W)  and  G^  G  {U')  ,£  —  0,1  then 


X  *+1  <  cdl+^  {k  +  2)!  VA:  >  0. 

^3  H  (W)  -  ■i  ''  ’ 


(4.20) 
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Proof.  Let  8t  =  8-  f  A(5  and  =  e  +  f  Ae,  where  AS  =  6  -  S'  and 

Ae  =  e'  —  e.  By  Ut  we  denote  Ue^^Sr  Then  U  =Uq  CU\  CUi  -  •  •  CUk  =  W . 

If  the  condition  (Bl)  holds,  the  application  of  (4.4)  leads  to 


V(Uh) 


+  E 


+  M{k)  ||u, 


<  Co  ^ll/r3llL^(i/^_j)  +  ^  lh^3  L^(Wk_i)) 

<ciMik)  (ii/iil^..)+,5,  II^1Ih-v.)) 


where  M  (k)  =  kMo  =  Armax  (^,  ;^).  The  argument  above  can  be  carried  out 
for  all  1  <  m  <  k+1.  Hence  e  (if) ,  1  <  m  <  A:  +  1,  and  by  the 

mathematical  induction  it  can  be  shown  that  for  0  <  s  <  A: 


(4.21)  < 


+  E  Gi 


+m)  f=0,l 


3  H  (U^-,+m) 


Then  (4.18)  follows  from  (4.21)  immediately. 

If  the  condition  (B2)  holds,  we  can  analogously  prove  by  mathematical  induction 
that  u^s+i  G  (if)  for  0  <  5  <  A:,  and 


(4.22) 


TT*/,,  \ 

^3  rl  (Uk-i-i) 


<  Co  {\\M\L,,,iUk-,)  +  ,5^ 


+  C^M-{k)  I  Co  ( \\f^s..-r. 


+  /xf+i-” 


+  E  bl 


II L/3j2  (^fc-5+m+l )  £=0,1 


TT2-f,2-£ 


ni:r'iUk-s 


+crv«  (i)  (ii/iil,„,«., + ,E,  I|g1Ih--V)  +  • 

We  shall  prove  (4.22)  for  =  0,  f  =  0, 1.  The  proof  for  the  case  that  0  is 
similar  to  what  follows.  (4.22)  holds  for  s  =  0  due  to  (4.5)  of  Theorem  4.1.  Suppose 
it  is  true  up  to  s,  then  applying  (4.5)  to  we  obtain 


<  Co|l/.|«  IIl^, +  C,M(k){C4S..*, 

+  _£  C„”M”(t)(Co||/,.+.--|lL,.,, 

+crM'«(t)(ll/llL,„,„.,  +  IN«llL=,„.,)} 

=  II IIl,,, (*)  {<^"ll/.|«-"llL,, 

+lk'S«-"  II^, (*')  (ll-f  IIl,„(«.)  +  11“*.  llL*(tt.))  • 


Hence  (4.22)  holds  for  (s  +  1),  and  then  we  complete  the  induction. 

If  /  €  B^(W')  and  then  there  are  some  da  >  1  and  Ca  >  1  such 

that  for  k>Q 


(4.23)  <  Czm, 

|ff'.  ^  (*  +  2  -  t)'- 

Substituting  (4.23)  into  (4.22)  and  noting  that  M{k)  =  kMo  we  obtain 

(4.24) 

il“.5«llH'(»)  S  C{CzCz^*\k  +  l)\ 

+  E  (Co  +  \)C^M^dl*'-’"(k  +  1  -  my.e'  + 

m=l 

By  using  Sterling  formula:  fc!  =  k^e~^\/2'Kk  ^1  +  0  (l))?  we  have  for  1  <  m  <  A: 

A:’”  (A:  +  1  -  m)!  <  CA:*"  (A:  +  1  -  m)*'+^"’”  e-('=+^-”‘y27r  (A:  +  1  -  m) 

(4.25)  <  C  (A:  +  1)*+^  e-(fc+i)y27r  (A:  +  l)-  e"* 

<  C(A:  +  l)!e™. 

We  have  by  substituting  (4.25)  into  (4.24) 
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TT* 


<  C I CoC3dl+\k  +  2)!  +  E  (Co  +  l){CoMoerds'^+'^-^{k  +  1)! 

m=l 


+(C'oMoe)''+i(fc  +  1)!} 

<Cdl+^{k +  2)1 


where  ds  =  max  (ds,  CoMoe^ . 

4.2  Regularity  of  high-order  derivatives  with  respect  to  the  direction 
perpendicular  to  the  edges 

We  now  turn  our  attention  to  the  regularity  of  high-order  derivates  with  respect 
the  variables  other  than  x^.  Let  y’l  (^3)  and  9^2  (^)  he  the  C°°  cut-off  functions 
defined  in  (4.13),  and  let  u  (x)  =  ifi  (xa)  (r)  u  (x)  where  u(x)  is  the  weak  solution 
of  (3.1).  Then  v  satisfies 


(4.26) 


—Av  =  f  mW  =  Ue'fi 

V  Ir=e'  =  ^L3=±(l-fi')  = 


v\e^o=~g^  =  G°l^^, 


where  =  (ft  (xa)  (p2  (r)  (xa)  if  2  (r)  and  f  —  fi  (xs)  V^2(r)/  -1-  h, 

h  =  2(Vi2U  •  Vi2(/72(r))v?i(x3)  -I-  ufi{x?)Ai2f2{r) 

+2u^M^^)^'ii'^)  +  ^^<^i(^3)</>2(r) 

with  V12  =  (^ — ,  ^12  =  ^”2)’  Obviously  /  €  L/j,2(ZY'),  &  € 

^dxi  dx2  dxj  dx^ 

(l/'),£  =  0, 1.  Furthermore  v,f  and  G^  vanish  for  r  >  e'  or  [xal  >  1  - 
and  with  the  constant  Mq  given  in  (4.7). 

(4.27a)  <  G  ^||/|Il,.2(^/')  +  ^0  \MViU')  +  X  1  vIL^w'))  ’ 


(4.27b) 


C°  tt2,2 


,  <c(e_m^-|(3«Ih>.,  +  IG”Ih“,„, 


<  CMq  *  ; 


(4.27c) 


<  CMoWG^Wjjii  . 


Iff  eV  (W)  and  G^  €  {U%  then  f  eV  (W) ,  G^  G  {W),  and 


(4.280.)  <  C  (||/IIl=,„,,  +  A/„  E  Mi 

\  8 — 0,1 


V(W)  I  ’ 


(4.28b) 


<  C  E 


Ye-\u') 


\}r(u') 


<  CMl 


i^{u')  =  [u£n\u')  I  u|,=e' =  wU=±(i-«')  =  o} 


ii]y{u')  =  \ueii\u')  I  wUo  =  o}- 

Then  n  -  G°  €  Hp  {W)  and  satisfies  the  following  variational  equations 

(4.29)  /  Vn  •  Vwdx  =  /  /dx  +  /  g^wdS,  Vi«  G  (W') . 

^  ’  Ju>  Jw  JT2f)dU' 

We  now  extend  u,  /  and  (5^  into  Q'^  x  IR^  by  zero  extension  outside  Id'.  Then  for 
almost  every  x  =  (xi,  X2)  G  Qe',  v  (x,  •)  G  (®'^)  ’  /  ’)  ^  ^  ^ 

(iR^)  ,£  =  0,1.  Let  .7^  denote  the  Fourier  transform,  namely,  for  admissable 
function  w 


1  '  \ 

wix,X)=J^{w)  =  -y=  u;(xi,X2,X3)e"*^®^dx3,  A  G  (-00, 00). 

yZTT  J-oo 

Then  v  =  (v) ,  f  =  [f)  ,  &  =  :F  (G^)  J  =  0,1  exist  and  v  solves  the  following 

problem 


(4.30) 


-Alt)  +  A^v  =  /,  in  Qe’  =  Q', 

u|r=e'  ~  0, 

A 

^1^=0  =  =  G°|0=O,  ■^|e=w  = 


Let  Hi,  {Q')  =  {w€H^  {Q')  1  w  =  0}.  Then  u  -  G°  G  Hi,  (Q')  and  satisfies 
the  variational  equation 


(4.31) 


■  V.tii  +  di=  I  fwdx  +  I  g'wds,yw  €  UpiQ') 

'  '  /  JQ'  J^2 


where  72  =  r2  D  dQ'. 
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We  shall  introduce  the  weighted  Sobolev  space  {Q')  defined  in  [].  For  integer 
jfc  and  fc  >  ^  >  0,  {Q')  is  the  completion  of  C°°-  function  in  the  norm 


(4.32) 


luWlrk 


mjQ') 


with  r  {x)  =  \x\  =  (xf  +  x^)^  ^  and 


/  r  (5) 


{-)  = 


for  |a'|  =  cii  +  oi2  > 
for  la'I  <  £. 


As  usual  we  shall  write  ((^0  =  iQ')- 

Lemma  4.2  Let  u  €  (il)  be  the  weak  solution  of  the  problem  (3.1)  with 
feLi3  (fi)  and  (11)  ,£  =  0, 1.  If  €  (0, 1)  satisfies 


(4.33) 


^12  >  1  ~  '^12)  ^^12 


Then  the  weak  solution  v  of  the  problem  (4.30)  belongs  to  {Q')  and 


(4.34)  <  C 


(4.35) 


<  C  ^1  +  |A|‘ 


1-2  +  Z)  4" 

fco.i  “/3i2  1 


L  m  I  ■ 


Proof.  Since  u  €  [Q'  x  1R^),  /  €  [Q'  x  IR^)  and  &  e  (Q'  x  IR^), 

V  €  ((5')?  /  €  and  e  H^72^’^“^ {Q')  for  almost  every  A  €  IRF 

Because  v  is  the  weak  solution  of  the  problem  (4.30),  by  Theorem  2.1  and  Remark 


of  [2],  V  €  Hj’"  (go  ,  and 


^  ^  ( 1|/ -  ^'"1IL«5')  +  Pfnxr'm 


with  C  independent  of  A.  This  leads  to  (4.34)  immediately. 

We  now  assume  that  =  0.  By  Lemma  2.10  and  2.11  of  [2]  we  get  for  any 

wen^  (go 


I  /  IIHIhOq') 


and 


Letting  w  =  v  and  substituting  the  inequalities  above  into  the  variational  equation 
(4.30)  we  obtain 


(4.36a)  + 

and 


/ 


+ 


L*,(«')  -  ir 


(4.36b) 

Here  we  used  the  inequality 

The  combination  of  (4.34)  and  (4.36)  leads  to 


161'" 


■k«,8.,  +  |A|''l|V.6|||,a„,,  +  |Ar||6||},.,,,, 


(4.37) 


<C 


(l  + 1^1")  (||/ 


L^12(<3') 


+ 


If  (7°  0  and  =  f  =  0,  setting  ?h  =  u  —  G”  and  applying  the  result  above  we  get 


which  together  with  (4.37)  implies  (4.35)  and  then  completes  the  lemma.  □ 

Theorem  4.3.  Let  u  €  H^(fl)  be  the  weak  solution  of  the  problem  (3.1)  with 
/  e  L^(H)  and  G‘  €  H^"'’"-'(0),  /  =  0,1.  In  addition,  f,,  G  L^,,(W')  and  G‘^^  € 
H^“J’^“*(Z7'),  /  =  0, 1  with  /3i2  satisfying  (4.33).  Then  u  G  H^’j^(Z7),  and 


(4.38) 


-  + ,5., 


^C'{||/||L,^^(n)  +  ll/-3llL,,,(ev') 


+  E  (Ill'll 

/=0,1 


H 


2-1, 2-1 
012 


(n) 


H 


2-1, 2-1 
^12 
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(4.39) 


S  c{ll/llL,„,«.,  +  ,5,  IG'Ih--'(«.,  +  ^^«II«’IIl»,«., 


where  Mq  is  a  constant  given  by  (4.7). 

Proof.  Since  /  €  and  G‘  €  Hj7j’^"'(770,  /  €  x  IR^)  and  G‘  € 

X  IRi),  /  =  0,1.  Hence  /  G  L^a2(<5')  and  G'  e  for  almost 

every  A  G  IR^,  and 


(4.40a) 

and  for  Z  =  0, 1 
(4.406) 


L^12(W')’ 


/oo  ^ 

.JIG'll  =  IIG'IIL 


'^12 


By  Lemma  4.2,  v  =  .F(v)  €  ^(^0^  and  (4.34)-(4.35)  hold.  If  /cm  G  L/3j2(7/') 

and  G^m  G  H^7J’^~'(77')  for  m  =  0, 1  and  Z  =  0, 1,  due  to  (4.34),  we  have 

by  (4.40) 
by  (4.27) 

This  leads  immediately  to  (4.39). 

Similarly,  due  to  (4.35),  we  have 

/OO 

(4-«)  <  cjy  +  |Ap)(ll/IIL,„„,,  +  E  IIG'llk--,,.,  +  |AniG”llb,,,,)<'A 
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where  =  <l)ii^3)(j>2{r)fx^  and  =  <l>i{x3)<l>2ir)G‘^^.  Therefore  we  have  estimates 
of  fxi  and  Gipj,  which  are  similar  to  those  in  (4.27)  except  that  /,  &,Ur  and  are 
replaced  by  fxi-,&x^,Urx3  and  u^>+i  .  Combining  these  with  (4.27)  and  (4.41),  we 
obtain  (4.38).  ^ 

Remark  4.1  Although  there  are  similarities  between  the  regularity  of  solutions  in 
the  edge  neighborhoods  and  those  of  solutions  for  the  problems  in  polygonal  domains 
in  IR^,  there  are  substantial  differences.  For  the  problems  on  a  polygonal  domain 
ft  in  IR^,  the  boundary  values  problem  of  Poisson  equation  realizes  an  isomorphism 
Hj+^’^(ft)  Hj’®(ft)  X  X  >  0,  but  it  is  not  true 

for  the  problem  on  a  polyhedral  domain  ft  in  IR^,  namely,  the  conditions  that  /  £ 
H^+2’^(ft),G'  €  =  0,1  are  not  sufficient  to  guarrantee  the  solution 

u  €  ° 

Theorem  4.4  Let  f{x)  G  and  ^  =  0, 1.  If  G  ) 

and  e  hJ;''/"'’^"'(W')  for  m  =  0, 1  and  /;  >  0  with  /?i2  satisfying  (4.33),  then 
the  problem  (3.1)  has  a  unique  solution  (weak)  u  €  H^(ft)  n  H^^^’^(W),  and  for 
\a\  ^  k  -\-2 


(4.42) 


+ll/llL,(a,+  nJ|G'llHr-(«)}- 


Furthermore,  if  /  €  and  G‘  €  =  0,1,  then  u  €  BlJU),  and 

there  are  some  constants  C  >  1  and  dj  >  1  such  that  for  all  a 


(4.43)  £  Cd-al 

Proof.  The  assertion  for  [aj  =  2  is  ture  owing  to  Theorem  4.3.  Let  6*  =  |(d  +  d') 
and  e*  =  |(^  +  S').  (4.42)  and  (4.43)  hold  over  U*  =  Ue*,s*  for  |a|  =  >  2  and 

cei  +  a2  <  2  due  to  Theorem  4.2.  It  remains  to  show  (4.42)  and  (4.43)  for  a  with 

.  e*  —  e 

ai  +  a2  >  2.  To  this  end,  set  Ui  =  We^,5i,0  <  «  <  fc,  with  c*  =  e  +  ^ ^ — ,  o*  = 

C  _  C* 

S  —  i — - — ,  so  that  7/o  =  S/c,s  ~  and  Lik  =  l^e*,s*  ~  SI* .  Let  v  —  ,  s  + 1  ^  fc, 

k 

s,t  >  0.  Then  v  satisfies 


(4.44) 


-An  =  /.,,  =  r-2(rVx‘)r^,  in 

^1^=0  =  r®G°s3,tlfl=o  =  Gjj|e=o, 

-^\e=ujn  =  r®"^(rGjt)r4|fcwi2  =  G\^\e=w32^ 
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Obviously  fs,t  e  Lp,,  (W),  Gi,,  G  1  =  0,1  and 


Due  to  Theorem  4.2,  we  have  for  t  =  k  and  s  =  0 


(4.45) 


+II/IIl„o,  +  ,5  IIG'llm— (»,)■ 


Suppose  (4.45)  holds  up  to  (s  —  1)  with  0  <  s  +  t  <  k.  Then  the  application  of 
(4.39)  of  Theorem  4.3  to  the  problem  (4.44)  gives  us 


E  11^^ 

|a|=|a'|=2 


t2 


L  («,_,) 


(4.46) 


+ll(''‘“rT;+”)'-llL’(«..,.,)  +  II II L’ 


By  the  assumption  of  the  induction,  we  have  for  m  =  0, 1 


<  C(||r-+'>-«,,....«  IIl",„,__^.,  +  Ik-”'’-..,..,*™  IIl=|„._..,|) 

(4-«)  <  C(||/..||hJ..  + ,  E_  l|Gi.llH”-''-(«,,  + 

+  E  l|G||„.-,,.-,,„j 


'Hr"  (H)' 


and  for  m  =  0, 1, 2 


<  ll-•-^■^'’■■>‘...r-llL’, 

(4-«)  <  C(||/.Jh«  +  ,5,  I|GL.IIh;«-'-,«'.  +  ll^«L,(a) 


+  E  IIGIIjj,. 

/=0,1  ^/3 


Combining  (4.46)-(4.48)  we  obtain  for  0  <  s  <  A:  and  s  +  t  <  k 


E  11$^; 

|a|=|«'|=2  ^  ^  < 


(4.49)  S  C'(||/.,|Ihw  +  II/IIlmo)  +  l|Gi.  IIh-’-'”' 


+I|G'|| 


Hr  ’  (0)' 
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This  completes  (4.42)  by  the  induction  for  a  with  |a|  =  A:  +  2  and  02  <  2. 

Now  it  remains  to  show  (4.42)  for  0.2  >  2.  For  a  with  |q!|  =  A;  +  2  and  0:2  >  2, 
let  w  =  r“*u^aiga2-22.“3.  Thcu 

-Au;  =  r“^“^(r^/e«2-2x“3)r“i . 

Noting  that 

Au;  =  r^^~'^'D^u  +  (2q:i  +  l)r°'‘“^Urai+iga2-2a;“3 

+ajr"*  ^^^01^02-2x33  +  ^°'*^r“l+2ea2-2j;“3  +  r"  ^U^aipa2-2x3  3+^' 


we  obtain 

<  ||r^^  ^(^^/0«2-"23,«3)r«i  +  (2<^  +  1)11^  ^  IIl^j2(^) 

(4-50)  IIl^i2(w)  +  lk"^^ir“i+30»2-2x“3 

+  ||r“>Vlfl«2-2x“3+HlL^^^(W)- 
Then  simple  induction  over  ot2  leads  to 

||r"i  ^Uroi5io2j,"3||L^^^(j^,j 

This  completes  (4.42). 

We  now  shall  prove  (4.43).  We  assume  that  =  0,  /  =  0, 1  for  simplicity,  and 
we  establish  the  following  estimates  by  mathematical  induction 

(4-51) 

5-f-2 

where  0  <  s  +  A  <  A,  Di  and  Dz  are  suitable  constants.  For  s  =  0,  (4.49)  holds  due 
to  (4.18)  of  Theorem  4.2.  Suppose  it  is  true  up  to  (s  —  1)  with  s  +  <  <  A.  Then  by 
application  of  (4.39)  of  Theorem  4.3  to  the  equation  (4.44),  we  have 
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2  2 

where  M,  =  max{^^,  — }  >  1.  By  the  hypothese  of  induction 


(4.53) 


<  C.{  E  E  ||/._i_,,,+„||r  .uMD^y- 

1=0  m=0  ^12  f 


(4.54a) 


'  ®3iiL  (Wfc_,+i) - 


<  ||r* 


(4.546) 


5-2 

+  E  Ih, 


m=0  ^3  ^  ) 

The  combination  of  (4.52)- (4.54)  leads  to 

lk®«r4+2®' |Il^^^(W)  <A  +  B-]rE, 


A  =  Co\\fsA\L,,,m  + 

+aE  E  E  ||/.-2-/,w'+™||r  + 

t'=0  /=0  m=0  ^12  ^ 

<  c.  E  E  I1/.-i..+~IIl,  iu'AD;’"k'-’", 

1=0  m=0  f 


B  =  C„C.{‘EJ|u„.+..|li,a|„.,C;C3”P-“+'M. 

+  E  E  i)!B3’“fc-’”+'-‘'M;-‘') 

t'=0m=0  ^^3  -L*  {C4  ) 


E  =(7„a{E 


V  — ^  II 

m=0  ^ 

2  5 


t'z=0m=0  3  ; 

5-f2 
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with  C*  >  Co  and  Di  >  4(7*  max(2  +  M*,  dz).  This  completes  the  induction. 

If  /  €  (?'  =  0,  /  =  0, 1  then  by  (4.19)  of  Theorem  4.2,  we  have  for  a  with 

|q;|  =  ai  +  02  <  2, 

(4.55a)  ^ 

and  we  may  assume  for  all  a  that 

(4.55(.)  CiZ)fD?D?a!. 

Substituting  (4.55)  into  (4.51),  we  obtain  for  s  +  t  =  k 

<  C4t  i  CrD{Dl{s-l)\{t  +  m)\k‘-^ 

l—O  m=0 

+  E  CiDl-^^Dlk^-^+\t  +  m)\+  *E  +  m)!} 

m=0  m=0 

<  a{4(7iT)fI>^t!/b*  +  2CiDl+^DltW+^  +  2CiDl+^  Dlt\k‘^+^} 

<  C2Dl+^Dik^+*+\ 

Note  that  by  Sterling’s  formula 

j^s+i+2  ^  )^k+2  <  2)!e*+y27r(lb  +  2)  <  Czik  +  2)!(2e)*+2 

which  implies 

l|r*t<,w..|lL,„(«.)  <  C,D‘*^DU2ef*'‘(k+2y.  <  C^iJ+VsiHs  +  2)!. 

Hence  (4.43)  is  proved  for  a  with  02  =  0,di  =  4eZli  and  dz  =  AeDz-  In  the 
same  manner  (4.51)  and  (4.56)  can  be  proved  for  a  with  02  <  2  except  that  v  = 
r"iu,.ai^a22;“3,  |q:|  <  k  and  02  <  2  satisfies  the  equation  (4.44)  and  the  estimates 
(4.45),  (4.51)  and  (4.56)  instead  of  u  =  with  s  1  <  k.  It  remains  to  prove 

(4.42)  for  a  with  02  >  2.  Suppose  that  (4.43)  holds  up  to  (02  —  1)5  then  by  (4.50) 
and  (4.55)  we  obtain 

<  +  C4{{2a  +  +  l)!(a2  -  2)103! 

+o2d“M^^-2d|®oi!(o2  -  2)103!  +  dr+"d“=-''d“^(oi  +  2)(o2  -  2)!o3! 

+d?>d^^-V^^+2oi!(ox  -  2)!(o3  +  2)!} 

<  C'4d"o! 

where  C4  >  Ci  and  dz  >  2ma.x{di,  D2,dz).  This  completes  the  induction.  Then 

(4.43)  holds  for  all  o,  and  u  € 
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Remark  4.2  If  f{x)  e  L/3(fi),  G  with  Pu  =  0,  /  =  0, 1,  then  (4.4) 

of  Theorem  4.1  and  (4.34)  of  Lemma  4.2  implies  that  u{x)  €  with 

satisfying  (4.33),  and 


Ih|^(«,  s  c’CII/iil,,.,  +  E  IIG' 


H  2-1, 2-r^^). 
3 


Further,  if  f{x)  €  Hj;°(i/')nL^(fi),  G'  €  H5+/-'-2-'(i/')nHj"'-^"'(a)  with  ^2  =  0, 
then  u{x)  G  with  ^12  satisfying  (4.33),  and 


H2— 1,2— I 

/3  (^) 


+I|G"| 


/3i2  (")■ 


□ 

As  a  consequence  of  Theorem  4.4  and  Theorem  2.1,  we  have  regularity  of  solution 
in  the  countably  normed  space  with  weighted  C*'—  norms. 

Theorem  4.5  If  f{x)  €  B^jj(W')  fl  Lp{f}),  G‘{x)  G  (W')  fi  H^‘”^’^~*(fl),  then 
the  weak  solution  u{x)  of  the  Poisson  equation  (3.1)  belongs  to  nH^(n).  □ 

Remark  4.3  In  special  case  that  the  data  functions  are  analytic  or  piecewise 
analytic,  namely, 

(i)  function  /  is  analytic  in  fl, 

(ii)  =  0, 1,  are  analytic  on  every  face  f,-  C  r°  and  f  j  C  f^. 

Then  the  solution  u  of  the  problem  (3.1)  belongs  to  and  H 

H\n).  ° 

Remark  4.4  The  regularity  described  by  the  countably  normed  space  with  weighted 
C*— norm  implies  the  pointwise  estimates  of  the  derivatives  of  solution  of  all  orders, 
namely,  for  x  G  1q;|  =  A:  >  0, 


(4.57)  \D°‘Hx)  -  u(0,0,a;3))|  < 


and 

^u(0,0,X3)  <C4kl 

In  many  applications,  for  instance,  the  error  analysis  of  the  p  and  h—p  version  of 
the  finite  or  boundary  element  method,  we  prefer  to  use  the  pointwise  estimates  of 
the  high  order  derivatives  of  solution  instead  of  the  weighted  Sobolev  norms  of  high 
order  derivatives.  By  using  estimates  (4.57)  and  (4.58)  we  have  shown  in  [8,9,16] 
that  the  approximation  to  functions  belonging  to  converges  exponentially 

by  properly  designed  piecewise  polynomial  spaces.  □ 


(4.58) 
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